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Abstract Owing to the Kubo relation, the shear viscosities of pionic and nu¬ 
cleonic components have been evaluated from their corresponding retarded cor¬ 
relators of viscous stress tensor in the static limit, which become non-divergent 
only for the non-zero thermal widths of the constituent particles. In the real-time 
thermal held theory, the pion and nucleon thermal widths have respectively been 
obtained from the pion self-energy for different meson, baryon loops and the nu¬ 
cleon self-energy for different pion-baryon loops. We have found a non-monotonic 
momentum distributions of pion and nucleon thermal widths, which have been 
integrated out by their respective Bose-enhanced and Pauli-blocked phase space 
factors during evaluation of their shear viscosities. The viscosity to entropy density 
ratio for this mixed gas of pion-nucleon system decreases and approaches toward 
its lower bound as the temperature and baryon chemical potential increase within 
the relevant domain of hadronic matter. 


1 Introduction 

A strongly interacting matter is expected to be produced, instead of a weakly in¬ 
teracting gas, at RHIC and LHC energies as the shear viscosity of the matter thus 
produced is exposed to be very small. This was concluded by the hydrodynami- 
cal simulations [ll[2l[3jlll[5l|6] as well as some transport calculations [71[8ll9l[T0] to 
explain the elliptic flow parameter observed at RHIC and LHC. According to the 
investigations conducted in Refs. the shear viscosity to en¬ 

tropy density ratio, 77 /s may reach a minimum in the vicinity of a phase transition, 
which is also indicated by some lattice QCD calculations [18[ll9ll20) . The minimum 
value of 77 /s may be very close to its quantum lower bound, commonly known as 
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the KSS bound |21| . Owing to these interesting issues, a growing interest in the 
microscopic calculation of shear viscosity for the QGP phase |221I231[2HI251I261I27| 
and hadronic phase pHl[^[3ni[5Tl[5^[55l[53l[55l[5Bll571[551[5^HDl[inH^[351fi31fi51HB] 
has been noticed in recent times, though the transport coefficient calculations of 
nuclear matter started somewhat earlier |471l481fi^l5f)ll511l52| . Importance of know¬ 
ing the explicit temperature dependence of shear viscosity for hadronic phase has 
been pointed out in a recent work by Niemi et al. [3]. They have shown that the 
extracted transverse momentum pt dependence of elliptic flow parameter, V 2 {pt), 
of RHIC data is highly sensitive to the temperature dependent 77 /s in hadronic 
matter and almost independent of the viscosity in the QGP phase. 

Being inspired by this, we have studied the shear viscosity of pionic medium [53] 
and then extended our study to the pion-nucleon system |54| . To calculate shear 
viscosity of the pionic component via Kubo relation |551I56| . the thermal corre¬ 
lator of viscous stress tensor for pionic constituents has to be derived, where a 
finite thermal width of pion should be included for getting a non-divergent value 
of correlator in the static limit |32ll401f^ . In Ref. [53], the pion thermal width 
is estimated from the pion self-energy for different mesonic loops, which are ob¬ 
tained in the formalism of real-time thermal field theory (RTF). Similarly, the 
thermal correlator of viscous stress tensor for nucleonic constituents is obtained in 
Ref. |5l| to calculate the shear viscosity of nucleonic component, where different 
pion-baryon loops are taken into account to determine the nucleon thermal width. 
Now, in the two component nucleon-pion system, pion propagation may also have 
some baryonic fluctuations besides the mesonic fluctuations. This contribution, 
which was absent in our previous studies |531I54| , is considered in the present work 
to revisit our shear viscosity calculation for two component nucleon-pion system. 

In the next section, the formalism of shear viscosity for pionic and nucleonic 
components are briefly described, where their corresponding thermal widths are 
discussed in three different subsections. In the subsection 12.11 we have elaborately 
deduced the pion thermal width in baryonic medium by calculating the pion self¬ 
energy for different baryonic loops in the formalism of RTF. In the next two sub¬ 
sections, the relevant expressions of pion thermal width due to different mesonic 
loops and nucleon thermal width from different pion-baryon loops are briefly ad¬ 
dressed as the detailed deduction of these expressions are already provided in the 
previous studies |531I54| . The numerical outcomes are discussed in Sec. |3| and in 
Sec.m we have summarized and concluded this article. 


2 Formalism 


Owing to the famous Kubo formula |55II56| . the spectral function of two point 
viscous-stress tensor, can determine the shear viscosity in momentum space 
by the standard relation |32| 

lim ^ /'d'‘a;e"‘''"^([7r,j(a;),7r"^(0)])/3 , (1) 

20 <7o,q^O go J 


where { 0)0 for any operator O denotes the equilibrium ensemble average; (0)0 = 


Tr 


eZfZO 

Tre-Z^-H ■ 
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The simplest one-loop expressions of Eq. m for pion and nucleon degrees of 
freedom are respectively given below m 

(wfc ){1-b nfc(ajfc)} (2) 

and 

VN = f [nt+ nl (wf ){1 - (wf)}] . 

foTT Jp 

(3) 

Their schematic diagrams are shown in Fig. [Ha) and Ha) respectively. Hence, 
adding the pionic and nucleonic components, we get the total shear viscosity 

VT = Vtt + VN ■ ( 4 ) 

In the above equations, nk(u;k) = — 1} is Bose-Einstein (BE) distribution 

of pion with energy = (k^ -|- whereas -|- 1} are 

Fermi-Dirac (ED) distributions of nucleon and anti-nucleon with energy uj^ = 
(k^ -|-m^)^/^. The corresponding thermal widths, and Tjv for pion and nucleon 
can be defined as 

B M 

= = a;fc,k)/m,r - (^o = wLk)/m^ 

B M 

(5) 

and 

■Tjv = = -^ImZ'^(^B)(fco = wf ,k) (6) 

B B 

respectively, where (fc) is pion self-energy for different nucleon-baryon 

(NB) loops (shown in Fig. [DJc) and (d)), (fc) is pion self-energy for dif¬ 

ferent pion-meson (ttM) loops (shown in Fig. [Ht')) and ^^(.,rB)(^) i® nucleon 
self-energy for different pion-baryon (ttH) loops (shown in Fig. Hb)). The super¬ 
script R stands for retarded component of self-energy and subscripts represent the 
external (outside the bracket) and internal (inside the bracket) particles for the 
corresponding self-energy graphs as shown in Fig. Hb), (c), (d) and Fig. Hb). 
Adoption of finite thermal widths TW and TW in Eq. 0 and © respectively is a 
very well established technique mM, which is generally used in Kubo approach 
to get a non-divergent value of the shear viscosity coefficient. In this respect this 
treatment is equivalent to quasi-particle approximation or relaxation time approx¬ 
imation. Again, this one-loop expression of 77,7 or r/N from Kubo approach |32II401 
I58| exactly coincides with the expression coming from the relaxation-time approx¬ 
imation of the kinetic theory approach EZIIIHIIIIIEH]. Hence, the thermal width 
of medium constituent plays a vital role in determining the numerical strength of 
shear viscosity of the medium. 

Next we discuss the calculations of thermal widths from different one-loop 
self-energy graphs as shown in Fig. o and ©• 


/37^ /■°° dkk® ^ 
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Fig. 1 The diagram (a) is a schematic one-loop representation of viscous-stress tensor for the 
medium with pionic constituents. The double dashed lines for the pion propagators indicate 
that they have some finite thermal width, which can be derived from the pion self-energy 
diagrams (b), (c) and (d). The diagram (b) represents pion self-energy for mesonic (ttM) 
loops. Direct and cross diagrams of pion self-energy for NB loops are represented by (c) and 
(d) respectively. 



Fig. 2 The diagram (a) is a schematic one-loop representation of viscous-stress tensor for the 
medium with nucleonic constituents and the diagram (b) represents nucleon self-energy for nB 
loops. 


2.1 Pion thermal width for different baryonic loops 

Let us first concentrate on the pion self-energy calculations for different possi¬ 
ble baryon loops i.e During propagation in the medium, pion propaga¬ 

tor can undergo different intermediate NB loops, where B = Z\(1232), ^"*(1440), 
N*(1520), 7V*(1535), zi*(1600), Z\*(1620), 7V*(1650), Zl*(1700), A^*(1700), N*(1710), 
N*(1720) are accounted in this work. The masses of all the 4-star baryon reso¬ 
nances (in MeV) are presented inside the brackets. The direct and cross diagrams 
of pion self-energy for NB loops have been represented in the diagrams [TJc) and 
(d). 

In real-time formalism of thermal field theory (RTF), self-energy becomes 2x2 
matrix with 11, 12, 21 and 22 components. From any of the components, one 
can found the retarded part of self-energy, which is directly related with physical 
quantity - thermal width (inverse of thermal relaxation time). Let us start with 
the 11-component of in-medium pion self-energy for NB loop : 


nl\NB){k) —i T 

a=-l, + l 


^L{k,l)EU{l)EU{l 


ak) 


(7) 





































Shear viscosity from mesonic and baryonic fluctuations 


5 


where and E]^{1 — ak) are scalar parts of the nucleon and baryon propa¬ 

gators respectively at finite temperature. In RTF this expression is as follows 

= T 2 - ~l , ■ - + nie{-lo))5(i^ - ml) , (8) 

— mj^ -H irj 

where nl{ijj^) = -)- 1} are the FD distributions of nucleon and 

anti-nucleon for energy = ( 1 ^ -|- and /tjv is the chemical potential 

of nucleon which is supposed to be equal with the chemical potentials of all the 
baryons considered here. The two values of a in Eq. © correspond to the direct 
and crossed diagrams, shown in Fig. [1] (c) and (d) respectively, which can be 
obtained from one another by changing the sign of external momentum k. 

Let us first discuss diagram (d) for which a = -t-1. Integrating Eq. ([7]) over 
and using the relation, 

lmn^(^jsfB){k) = ta.nh{/3ko/2)lmnlljsfB){k) , ( 9 ) 

the retarded component of the in-medium self energy (imaginary part) can be 
expressed as 

Im<^^)(*) = / {LrAul’uS 

^i[{l - nliojl) - nZ{uJu)}S{ko - ujj^ - uju) 

+ {n'l{uj^)-nt{u}^)}5{ko-uj^ + uj^)] 

-h L 2 [{-n;“(wi^) + n“(wf )}(5(fco - ^u) 

+{—l-t-n.; [uf') + nl{u^)}&{koLof' +w^) , (10) 

where nl(uj^) = g^j.g pp) distribution functions for baryon 

and anti-baryon with w® = {(1 — k)^ -|- and Li ^2 denote the values of 

L{Iq, 1 , k) for Iq = and —oj^ respectively. The different 5 functions in Eq. (1101) 
create the regions of different branch cuts in fco-axis viz. —oo to — {k^ -|- (mjv -h 
for unitary cut in negative fco-axis, —{k^ -|- (ms — mjv)^}^^^ to {k^ -|- 
{mB — for Landau cut and {k^ -|- (mjv -h to oo for unitary cut 

in positive fco-axis. In these different kinematic regions, the imaginary part of the 
pion self-energy becomes non-zero. Among the four terms in the right hand side 
of Eq. (11011 . the third term contributes in pion thermal width for baryonic loops, 
E-k(nb) because the pion pole (fco = w((,k) is situated within the Landau cut (0 
to {k^ -|- (ms — ) in the positive fco-axis. From the Eq. dS]), using the 

relation, 

TVcjvb) = -Inil7,f(^£j)(fco = Wfc,k)/m,r (H) 

and adding the relevant Landau cut contributions of both diagrams (c) and (d), 
the total thermal width of pion for any NB loop is given by 

N 

r,(^B)(k,T,MJv) = xL(Zo = -t^r,l,fco=wLk) [{-n+{ul) 

+71^ {u!u = )} + { — ril (wf^) -h (oj® = + Lof )}] , 

(12) 
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where 


with 


N 

Wi± 


rr2 

2m% 


{-ul ± |k| W^(mb)) , 


S‘^(NB) — ~ 


and 

W^7r(ArB) = (l — ^ . 


(13) 


(14) 

(15) 


Lagrangian densities of spin Jg = 1/2 and 3/2 baryons can be written as 


p _ X '' ^free . ^int 

'^Baryon — 2_^ |_'^S(Jb = 1/2,3/2) ^B{Jb = 1 !‘2,^ 12) 

B 


(16) 


where free parts of Lagrangian densities for baryonic fields with spin Jb = 1/2 
and Jb = 3/2 are 


C 

C 


free 

B(Jb = 1/2) 


free 

B{Jb=Z/2) 


B{Jb = 1/2) 


E f / / a o/3 - olB n I 1/ 

a-imBO- )1PB 


B(Jb=3/2) 
with cr“^ = ^ 


CK 8 

7 ,7 


(17) 


and their interaction parts are |59| , 

^mt ^ I I iPnO^^tt + h.c. for jf = I , ( 18 ) 

BT-n- L J ^ 

Here P stands for parity quantum numbers of the baryons. The coupling constants 
ttNB interactions are fixed from the experimental decay widths of B ^ Nn chan¬ 
nels [60]. They are f/m^ = 15.7, 2.5, 11.6, 1.14, 3.4, 1.22, 1.14, 9.5, 2.8, 0.35 and 
1.18 for Z\(1232), Ar*(1440), 1V*(1520), 7V*(1535), Z\*(1600), Z\*(1620), 1V*(1650), 
2\*(1700), 7V*(1700), 1V*(1710), 7V*(1720). Using (flSl) . we have found the vertex 
factors [S3: 


L{k,l) = -4 



8 

3to^ 


[2(fc ■ l)^ — a{k ■ l)k^ — {i^ -h mivms)], for Jb = 

7-^ [mjvTTiB + — a{k ■!)][(! • k — ak^)^ — fc^ms], 

TTItp J 

f jP 

forJs = - . 


1± 

2 


( 19 ) 
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2.2 Pion thermal width for different mesonic loops 


To calcnlate the mesonic loop contribntion of pionic thermal width we 

have evaluated pion self-energy for ttM loops, where M stands for a and p mesons. 
This contribution estimated in our previous work elaborately. Following that [53] , 

the expression of pion thermal width from the mesonic loops is given below 

r„{iTM)0<,T) = I - [ dujf L{lo =-ujf ,l,ko = 

167r|k|m^ 

- nu{u}^ =^k + wD} , (20) 


where ni, Uu are BE distribution functions of tt, M mesons respectively and the 
limits of integration are 


wz± = 


^7r(7rM) 

2771^ 


{-ul ± |k| IT^(^m)) , 


with 

and 


T(7rM) 

4 /o4 \l/2 


^^7r(7rM) (f /^'7T(7rM)'} 

Lagrangian density of pion, sigma and rho mesons can be written as 


jC = C 


free 


^free . ^free , ^int . ^int 

^<7 ' “T ^'K'KO “r •^TTTTO 


( 21 ) 

( 22 ) 

(23) 

(24) 


where free parts of Lagrangian densities for pseudo-scalar vr, scalar a and vector 
fields are are 

yc»free 2 2'\ 

- T^aO- ) 

4“ = - ml{p,p^)}, p^- = (a>^ - d-^pn (25) 

and their interaction parts are mMm 

= gppf,-'^ X d^-K 

= ^rriaTT ■ TTcr . (26) 

The coupling constant Pp = 6 and go- = 5.82 are fixed from experimental decay 
width [53] and physical masses of pion, sigma and rho mesons are taken as tUtt = 
0.140 GeV, TTia = 0.390 GeV and rUp = 0.770 GeV. Using (1261) we have obtained 
the vertex factors: 

2 2 

L{k,l) = -^^^, forM = cr, 

2 

= —-4- [fc^ (k^ — rrti') + — rrvi') — 2{ik ■ 1) m?p -|- l^}], for M = p . 

rrip \ c/ \ c/ 

( 27 ) 
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2.3 Nucleon thermal width 

In order to calculate the nucleonic thermal width rjv( 7 rB)j we have evaluated nu¬ 
cleon self-energy for different possible ttB loops, where B stands for all the baryons 
as taken in pion self-energy for baryonic loops. This contribution is rigorously ad¬ 
dressed in our previous work [SI] . Hence taking the relevant expression of nucleon 
thermal width for any ttB loop from the Ref. jSl] , we have 

rN{iTB){KT,HN) = J du}[L (^lo =-oji ,\,ko = uJk 

+ nu{cOu = <^k + i^D} I (28) 

where ni and Uu are BE and FD distribution functions for tt and B respectively. 
The relevant limits of integration in Eq. (1281) are: 

(-cuf ± |k| ITjv(.b)) , (29) 

with 

= m% -m% +ml (30) 

and 

fk"iv(7rB) = (l — ^ . (31) 

The vertex factors |54| : 

L{k, 1) = - I ^0 “ PmlmB | for jf = i , 

= -™'^s|(fco-/o + EmB)for = 5 

(32) 

can be deduced by using the ttNB interaction Lagrangian densities from Eq. m- 


3 Results and Discussion 

The detailed Landau cut contributions of pion self-energy for mesonic loops and 
nucleon self-energy for different ttB loops are investigated in the earlier Refs. |53| 
and m, where their corresponding contributions in the shear viscosity are also 
addressed. Now in the two component pion-nucleon system, another contribution 
to pion thermal width can arise from the pion self-energy with baryonic loops, 
which was not considered in our previous studies of the shear viscosity |531I54| . The 
main purpose of the present work is to include these baryonic loop contributions 
in the pion thermal width and to revisit the shear viscosity results. 

Let us first zoom in our attention on the T)r(jvB)- Fig- © represents the Landau 
cut contributions of different NB loops on the invariant mass axis Mfc, which can 
be numerically generated by replacing uj]^ = (k^ -|- in Eq. (1121) . For a 

fixed set of parameters k, T and /ijv, the r^{NB){Mk) for baryons B = Z\(1232), 
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Fig. 3 for different NB loops in their Landau regions, which contain the pion 

pole Mfc = m-TT, denoted by straight dotted line. 



Fig. 4 The on-shell thermal widths (upper panel) and mean free paths (lower panel) of pion 
for ttM loops, NB loops and their total are represented by dotted, dashed and solid lines 
respectively while dash-dotted line denotes the same results for nucleon component with all 
possible nB loops. 


Z\*(1650), A^*(1440), 7V*(1700) in the upper panel, B = 7V*(1520), Zi*(1620), 
A^*(1535), Z\*(1600) in the middle panel and B = A^*(1720), Z\*(1700) in the 
lower panel are individually presented in the Fig. ©■ The Landau cut regions, 
where r^{NB){Mk) for different NB loops have attained their non-zero values, 
are clearly observed in the invariant mass axis. As an example, for NA loop the 
Landau region is from M = 0 to (m^ — mjv) = 0.292 GeV. The straight dotted 
line denotes position of pion pole (i.e. Mk = TOtt), which indicates the on-shell 
contribution of r^(NB) for different baryonic loops. Here we identify NA loop as 
a leading candidate to contribute in the pion thermal width among all the baryonic 
loops. 
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Fig. 5 Same as Fig. along the /xjy axis. 



Fig. 6 Same as Fig. l|4]l along the k axis 


Adding the on-shell contribution of all NB loops, we get the total thermal 
width of pion for baryonic loops, which is plotted against temperature by dashed 
line in the upper panel of Fig. ©• Using Eqs. (1201) and (12811 , we have also generated 
the numerical values of total r^[TrM){T) (dotted line) and Ujv( 7 rB){T) (dashed line) 
by adding their corresponding loop contributions for a same set of input param¬ 
eters (k = 0.3 GeV, /rjv = 0). Following Eq. ([5]), the solid line represents the T 
dependence of total thermal width of pion, /^(T) after adding the mesonic and 
baryonic loop contributions. All of the U’s are monotonically increasing function 
but with different rate of increment. The corresponding results of mean free path, 
defined by A = k/(ci;fcU), are presented in the lower panel of the Fig. (|3]). Being 
inversely proportional to the thermal width, the mean free paths for all of the 
components monotonically decrease with T and exhibit divergent nature at low T. 
Along the /xjv axis, A’s (F’s) for all of the components also decrease (increase) with 
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Fig. 7 Temperature dependence of shear viscosities (upper panel) and entropy densities (mid¬ 
dle panel) for pionic (dotted line), nucleonic (dashed line) components and their total at two 
different nucleon chemical potentials: /rjy = 0 (solid line) and /rjy = 0.5 GeV (dash-dotted 
line). In the lower panel, the ratios of total viscosity to entropy density are represented as a 
function of T at same set of pjv’s, taken in the upper and middle panels. 


different rates as shown in the lower (upper) panel of Fig. ([5]). Here we see that in¬ 
dependent nature of pion thermal width (F.n-( 7 rM)) or mean free path for 

mesonic loops is transformed to an increasing or decreasing nature when the bary¬ 
onic loop contribution is added. Moreover, the divergence problem of \-k(nb){p-n) 
at low /xjv is also cured in the total mean free path for pionic component A 7 r(/ijv). 
A mild /t]v dependence of the nucleonic component is observed. 

At fixed values of T and /rjv, the momentum distribution of thermal widths 
(upper panel) and mean free paths (lower panel) for all of the components have 
been displayed in Fig. Being equivalent to the momentum distribution for the 
imaginary part of optical potential (see e.g. [MiEa), thermal width of pion for 
any mesonic or baryonic loop exhibits a non-monotonic distribution with a peak 
structure along the k axis. The mathematical reason can roughly be understood 
from the relevant Eqs. (Uni) and (USD as described in Ref. [60]. After adding the 
different NB loop contributions, each of which has similar kind of momentum dis¬ 
tribution with different numerical strength, we get a multi-peak complex structure 
of Hn-(iVB)(k). When we add it with the F^(.n.M)(k), which contains a dominating 
profile with one peak (due to np loop mainly), then a well behaving momentum 
distribution with less complex structures (solid line) is obtained. The FAr(k) (dash- 
dotted line) approximately appears constant with a mild reduction with k. Though 
we notice a divergent nature of \tt(nb)0^) out side the range of k = 0.1 — 1 GeV 
but the total A 77 in the entire momentum range remains non-divergent or finite 
with an well-behaved distribution. 

Using the total thermal width for pionic component, iV(k, T,/tat) in Eq. ([^ 
and for nucleonic component, Fjv(k, T,/tat) in Eq. (|3]), we have obtained shear 
viscosities ry 7 r(T,/ tat) and rj^iT, p^)- They are plotted by dotted and dashed lines 
respectively as functions of T and /tat in the upper panels of Fig. (|7|) and (| 8 |). After 
exhibiting a soft peak structure in the low T(< 0.1 GeV), rjTv{T) monotonically 
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Fig. 8 The dependence of same quantities as Fig. 0 at two different temperatures: 
T = 0.12 GeV (solid line) and T = 0.15 GeV (dash-dotted). 
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Fig. 9 Dashed and dotten lines show r]-,^ without and with baryonic fluctuations in pion 
propagation respectively. 


increases with a very mild rate in the high T(> 0.1 GeV) domain. When a mono- 
tonically increasing function r]]^{T) is added with this pionic component, the total 
shear viscosity rjT in high T domain enhances with slightly larger rate of incre¬ 
ment, which can be noticed by solid line in the upper panel of Fig. ©• Another 
curve of r]T{T) at = 0.5 GeV is shown by dash-dotted line, which faces a rapid 
increment after T = 0.06 GeV. The reason of this drastic enhancement can be well 
understood from the /rjv dependence of the two components 7777 (/r at) and 77jv{PAr)- 
The upper panel of Fig. m exposes a rapidly increasing function 77jv(/rjv) and a 
soft decreasing function r]Tr{fiN)- Remembering the Fig. ([5]), we can understand 
that the origin of soft decreasing nature of r]Tr{iJ-N) is coming from the baryonic 
loop contribution of pion as its mesonic loop contribution is independent of pjv- 
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Right panel of Fig. (|9]) is zooming this fact more distinctly, where we see how 
inclusion of NB loops in pion self-energy makes rj-^ deviate from its independent 
(dashed line) to dependent (dotted line) nature with . This is the main and dra¬ 
matically important contributions of the present article as an extension of earlier 
works |531I54| . After including it, providing a complete picture of shear viscosity 
calculation for pion-nucleon system is the main aim of this present investigation. 

The dotted lines in the left and right panel of Fig. m are exactly same as 
dotted lines in the upper panels of Fig. © and dSl) respectively. Still those curves 
are repeated for elaborating the effect of baryonic fluctuations in pion self-energy. 
As the phase space factor of Eq. does not depend on the fiM, so only thermal 
width T7r(/riv) controls on the pjv dependence of r/n- Now between two components 
r^i ttm) S'lid T)r(ArB) of , the latter one has only the dependency of hn as exposed 
in Fig. ®. The decreasing nature of is solely governed by the increasing 

(decreasing) nature of function T’ 77 (/rjv) Whereas, in case of Eq. ([3]), 

nucleonic phase space factor depend on /rjv so strongly that it makes r7jv(/rjv) be 
an increasing function after dominating over the opposite action of Tjv(/rjv) or 
Ajv(/rjv) on the 77jv(mjv)- 

Middle panels of Fig. Q and ([S]) represent the T and /xjv dependence of entropy 
densities for pionic and nucleonic components by following their ideal expressions: 


and 

SM = 4.p j + ^-fXN^n+{u:^) . (34) 

Using these numerical results of entropy densities (middle panel) as well as for 
shear viscosities (upper panel) for pionic, nucleonic components and their total, we 
have presented their corresponding ratios in the lower panels of the graphs, where 
straight horizontal (red) lines stand for KSS bound of the ratio. The decreasing 
nature of ratio is sustained for both /rjv = 0 (solid line) and /tat = 0.5 GeV 
(dash-dotted line) in the entire T axis. The former is dominating over the later 
in magnitude for T <0.12 GeV and then an opposite trend is followed beyond 
T = 0.12 GeV. Therefore, the ratio in the piv axis at T = 0.15 GeV (dash-dotted 
line) and T = 0.12 GeV (solid line) are exhibiting a nature opposite to each other 
(up to /rjv ~ 0.4 GeV), which can be observed in the lower panel of Fig. ([8]). 
Nevertheless, both of them increase in high /tat domain {/in > 0.4 GeV). Most of 
the earlier work msMm showed a reducing nature of ratio along the hn axis, 
which is also found in the present work up to T « 0.12 GeV but beyond T = 0.120 
it is not found. It indicates that our approach has some deficiency with respect 
to the earlier work |331l441fi5| . This deficiency may be the mixing effect of two 
component system |33| . which have been taken care in our further investigations 
and discussed in next paragraph. 

We have adopted a rough mixing effect which can generally be expected 
between two components of a mixed gas [31163]. From the Eqs. © and m , one 
can clearly notice that the phase space factors of rj-jr and ijn do not face any 
mixing effect of pion density, p.^ = 3 f and nucleon density, pat = 

4 J ( 2 Tr)^ )• Although their thermal widths F-^ and F^ contain this mixing 
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Fig. 10 Upper panel: Temperature dependence of shear viscosities of pionic (dotted line), 
nucleonic (dashed line) components and their total (at three different values of in presence 
of mixing effect. Lower panel: The ratios of total viscosity to entropy density vs T at = 0 
(solid line), 0.3 GeV (dash-double-dotted line) and 0.5 GeV (dash-dotted line). 



(GeV) 


Fig. 11 upper and middle panels show the dependence of same quantities as in Fig. (nnii 
at two different temperatures: T = 0.12 GeV (solid line) and T = 0.15 GeV (dash-dotted). 
Lower panel: The different points of (T, //at), where is approximately equal to the KSS 

bound. 


effect as they depend on thermal distribution functions of both, pion and nucleon. 
Following the approximated relation [3311541163] 

mix mix . mix /'otr\ 

Vtot = ?77r + ?7jv , (35) 




mix 

TV 


_ l]n _ 

i + (^) 


with 


( 36 ) 
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and 


mix 

Vn 


r]N 


1 + 


\PN J \ CTJVJV ) 


1+mjv / ’tT'tt 

2 


(37) 


where the cross sections of all kind of scattering are simply taken as constant with 
same order of magnitnde (i.e. ~ cp-kN ~ o’jvjv)- In presence of this mixing 

scenario, the T dependence of (dotted line), (dashed line) and their total 
77 ™))^ at p,jv = 0 (solid line), /tjv = 0.3 GeV (dash-dotted line) and /tjv = 0.5 GeV 
are shown in the npper panel of Fig. (nni. The corresponding results along the fx n 
axis for two different temperatures are presented in the upper panel of Fig. (HID . 
Lower panels of Fig. (nni) and middle panel of Fig. (ED are displaying the viscosity 
to entropy density ratios as a function of T (at three different values of /rjv) and /rjv 
(at two different values of T). One should comparatively notice the dash-dotted 
line in the lower panel of Fig. ([SD and the middle panel of Fig. ED, which are 
exhibiting an opposite nature in the low /ijv region. Hence the approximated T- 
/ijv range, where viscosity to entropy density ratio reduces, is transformed from 
(T = 0 - 0.12 GeV, /ijv = 0 - 0.5 GeV) to (T = 0 - 0.15 GeV, ^J,N = Q- 0.5 GeV) 
due to the mixing effect. Being closer to the earlier results mMm, specially 
the result of Itakura et al. [33], the mixing effect appears to be very important. 
Though the ratios in both cases, without and with mixing effect increase beyond 
the /rjv ~ 0.5 GeV but the conclusion of our results, based on the effective hadronic 
Lagrangian, should be concentrated within regions of 0.100 GeV < T < 0.160 GeV 
and 0 < /ijv < 0.500 GeV. The lower panel of Fig. ()11D shows the T-/rjv points 
where the ratios are approximately equal to its KSS bound. From this plot we 
can get a rough idea of T-pLjq region, where our shear viscosity calculations for 
the pion-nucleon system may be considered as a reliable estimation by using the 
effective hadronic model. 


4 Summary and Conclusion 

The present work is an extension of our previous studies |531I54| of the shear viscos¬ 
ity calculations for pionic |53| and nucleonic |5l] components, where pion thermal 
width due to different mesonic fluctuations and the nucleon thermal width due to 
different pion-baryon fluctuations are respectively considered. However, in the two 
component pion-nucleon system, pion thermal width may also be originated from 
different baryonic loops, which is not taken in our previous investigations iMEa. 
Gonsidering this baryonic loop contribution, we have addressed a complete pic¬ 
ture of pion and nucleon propagation via all possible meson and baryon quantum 
fluctuations at finite temperature and density, from where their corresponding 
contributions to the shear viscosities have been found. 

Following the traditional technique of Kubo relation |321l4nil58ll57| . the shear 
viscosities of pion and nucleon components can be deduced from their corre¬ 
sponding correlators of viscous stress tensor in the static limit, which will be 
non-divergent when a finite thermal width will be introduced in their free propa¬ 
gators. These finite values of pion and nucleon thermal widths have been estimated 
from the RTF calculations of the pion self-energy for different mesonic and bary¬ 
onic loops and the nucleon self-energy for different pion-baryon loops. Thermal 
width and its inverse quantity, mean free path for each component is numerically 
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generated as a function of the momentum k of the constituent and the medium pa¬ 
rameters, T and /ijv. They show very non-trivial momentum distributions, which 
have been integrated out by the Bose-enhanced and Pauli-blocked phase space 
factors of pion and nucleon, respectively, to calculate their corresponding shear 
viscosities. We have plotted the shear viscosity of each component and their total 
as a function of T and hn, where one can observe a distinct and important effect 
of pion thermal width for baryonic loops, which is the main finding of the present 
investigation to demonstrate a complete picture of shear viscosity calculation for 
pion-nucleon system. Actually the /ujv dependence is entering in the shear viscosity 
of pionic component via this baryon loop contribution of pion thermal width. This 
additional contribution makes the shear viscosity of pionic component reduce with 
/ijv and increase with T. By adopting a rough mixing effect of pion and nucleon 
densities between two components, we have tried to present a numerical estimation 
of total shear viscosity for a mixed gas of pion-nucleon constituents. Normalizing 
by the ideal expressions of entropy densities for pion and nucleon gas, we have ob¬ 
tained the viscosity to entropy density ratios for each component and their total. 
In the relevant region of hadronic domain, this ratio for the pion-nucleon 

gas mixture reduces and approaches toward its KSS bound as T or /ijv increases. 
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